Coulomb energy difference as a probe of isospin-symmetry breaking in the upper /p-shell nuclei 
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The anomaly in Coulomb energy differences (CED) between the isospin T = I states in the odd-odd N = Z 
nucleus ^"Br and the analogue states in its even-even partner ^'^Se has remained a puzzle. This is a direct 
manifestation of isospin-symmetry breaking in effective nuclear interactions. Here, we perform large-scale 
shell-model calculations for nuclei with A = 66 — 78 using the new filter diagonalization method based on the 
Sakurai-Sugiura algorithm. The calculations reproduce well the experimental CED. The observed negative CED 
for A = 70 are accounted for by the cross-shell neutron excitations from the /p-shell to the gg intruder orbit 
with the enhanced electromagnetic spin-orbit contribution at this special nucleon number. 

PACS numbers: 21.10.Sf, 21.30.Fe, 21.60.Cs, 27.50.+e 
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Isospin is a fundamental concept in nuclear and particle 
physics. The isospin symmetry was introduced under the as- 
sumption of charge independence of the nuclear force 
Historically, the study of this symmetry led directly to the 
discovery and understanding of quarks. However, it is well 
known that this symmetry is only approximate because of the 
existence of the Coulomb interaction and isospin-breaking in- 
teractions among nucleons, leading to small differences, for 
example, in the binding energy of mirror-pair nuclei and in 
the excitation energy of the same spin, 7, between isobaric 
analogue states (IAS) of the same isospin, T . 

A nucleus is a quantum many-body system with finite size, 
which generally shows two unique features in structure: the 
shell effect with the presence of strong spin-orbit interaction 
(I2I1 and the nuclear deformation associated with collective mo- 
tion f^. To properly describe these aspects in the frame- 
work of nuclear shell models, effective interactions must be 
involved. Thus the effects of isospin-symmetry breaking can 
manifest themselves through structure changes in the vicinity 
of the N line, providing information on the ^--dependence 
of the effective interactions. The effects have been extensively 
studied for nuclei in the upper sd- and the lower /;?-shell re- 
gions (see Ref. for review), where a remarkable agree- 
ment between experimental mirror energy differences (MED) 
and shell-model calculations has been found, allowing a clear 
identification of the origin of isospin-symmetry breaking in 
effective nuclear interactions. 

The Coulomb energy difference (CED), defined by 

CED(7) = r = 1 , r, = 0) - e,{j, r - 1 , r, - 1), (i) 

is often regarded as a measure of isospin-symmetry breaking 
in effective nuclear interactions which include the Coulomb 
force |5, 6]. In Eq. Ex{J, T, T^) are the excitation energies 
of IAS with spin J and isospin T, distinguished by different TJ. 
(the projection of T). 

More complex shell structures are expected for heavier 
mass regions. In the upper /p-shell, for example, abrupt struc- 




FIG. 1: (Color online) Experimental CED between the isospin T = I 
states in odd-odd N = Z nuclei and the IAS in even-even nuclei for 
mass numbers A = 42 — 78. Data were taken from Refs. 0-@,[Tl|]. 



ture changes along the N = Z line 10] and the phenomenon of 
shape coexistence are known. Being pushed down to 

the lower shell by the spin-orbit interaction |2], the §9/2 in- 
truder orbit and its interplay with the /p-shell orbits play a 
key role in the overall structure. However, the influence of 
these structure changes on CED has not been explored. In Fig. 
[T] we show the experimental CED for nuclei with mass num- 
bers A = 42 to 78, together with the very recent data of the 4j^ 
and 6^ states in A = 66 ifTlll . One observes that, while most 
of the CED values are positive, only the mass number A = 70 
has negative CED fd,L^. There have been suggestions to ex- 
plain such an anomalous behavior but without much success. 
In Ref. fl, it was suggested that this behavior could be at- 
tributed to the Thomas-Ehrman effect 1 13",^ of the loosely 
bound proton in ^"Br. However, as pointed out in Ref. fst], 
this suggestion cannot explain the observed positive CED in 
other systems such as A = 66, 74, 78 because these systems 
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have similar binding energy differences as in A = 70. In Ref. 
Jst], it was suggested that the anomaly may be qualitatively ac- 
counted for by small variations in the Coulomb energy due to 
shape differences in the pair of nuclei. The authors in Ref. |5] 
tried to argue that such deformation changes occur only for 
a prolate shape, and not for transition from an oblate ground 
state to excited configurations |5]. However, recent measure- 
ments and theoretical calculations indicate that both ^''Br and 
™Se are associated with an oblate shape 0121 llSll . Most re- 
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cently, the anomalous behavior has been suggested as due to 
different mixing of competing shapes within the members of 
the isobaric multiple! jl ill . 

In the past few years, experimental data for mirror nuclei 
above the doubly magic ^^Ni have become available. The 
MED in the A ~ 60 mass region were discussed |i , 16i,ll7ll . It 
was reported that the contribution to MED from the electro- 
magnetic spin-orbit term £/, is significant, while the monopole 
Coulomb term £// is not. Here, detailed variations in nuclear 
shell structure enter into the discussion. Energy shifts due 
to the Els term increase the gap between the 753/2, fsft 
the gg^2 orbit for neutrons but reduce that for protons. As 
a consequence, excitations involving these orbits contribute 
differently to MED. We have recently investigated MED be- 
tween the isobaric analogue states in the mirror-pair nuclei 



^^Se and ^^As using large-scale shell-model calculations 1 18|. 



The calculations reproduced well the experimental MED 
and suggested that the e/^ term provides the main contribu 
tion to the observed MED values. The negative MED values 
observed in the A = 67 mirror nuclei were explained by the 
single-particle-energy shifts due to the 6/^ contribution. We 
may thus expect that the Eis term plays a similarly important 
role for the anomaly in CED in the A = 70 IAS. 

To study CED and probe isospin-symmetry breaking in fi- 
nite nuclei, calculations for energy levels with an accuracy of a 
few keV are required. The current state-of-the-art shell-model 
calculations with the m-scheme using the Lanczos diagonal- 
ization method |20] can handle a matrix dimension of the or- 
der of 10^". However, it is difficult for the Lanczos method 
to obtain solutions for T ^ I states in odd-odd N = Z nuclei 
under the isospin-symmetry breaking with the Coulomb inter- 
action. As for instance, it is impossible to obtain the T ^ I, 
J — A state because of the presence of numerous T states 
below it. To solve this problem, a breakthrough in shell-model 
calculations is needed. Quite recently, the filter diagonaliza- 
tion method for large-scale shell-model calculations has been 
proposed by some of us |21], which is based on a new algo- 
rithm for diagonalization suggested by Sakurai and Sugiura 
1 2211 (called the SS method). The filter diagonalization proce- 
dure can overcome the difficulty mentioned above ||23ll . 

In this Letter, shell-model calculations are performed by us- 
ing the filter diagonalization method in the p/5/2g9/2 model 
space. We employ the recently-proposed JUN45 interaction 
1 2411 . a realistic effective interaction based on the Bonn-C po- 
tential and adjusted to the experimental data of the A = 63 ~ 
96 mass region. We try to solve a general shell-model eigen- 



H\^t)^ek\^,,), 
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where H is the shell-model Hamiltonian, and ek and jO^.) are 
eigenvalues and eigenfunctions, respectively. The spirit of the 
filter diagonalization method |21] is that in order to reduce 
the large-scale eigenvalue problem to a smaller one, one in- 
troduces moments pLp{p = 0, 1,2, • • • ) defined by the Cauchy's 
integral 



(3) 



where jy/) and |0) are arbitrary wave functions, e is a given 
energy for a target region, and Y is an integration contour 
Following the SS method li22,l , one can extract the eigenval- 
ues ei^{k E r) inside the closed curve F from these moments. 
In this way, the large-scale eigenvalue problem inside F is re- 
duced to 



Mx = XNx, 



(4) 



where M and are n x n Hankel matrices expressed by the 
moments /Xp and the eigenvalues Ao, ' ,Xn-i in the interior 
of F. After diagonalizing (|4]l, one obtains the eigenvalues 
eit — e + Xii inside the integration contour F (Details were dis- 
cussed in Ref. 12111 ). For the present case, we apply filter di- 



12D 7 agonalizations by taking the T \ Xj/J^ ' ) state as | Xj/) and | ^ ) in 
Eq. (|3]l for the odd-odd N — Z nucleus, where T is the lower- 
ing operator of isospin and | Wj^^) is an isovector T = 1 state 
with spin J in the even-even N = Z + 2 nucleus. Therefore 
r_|v/J=') is a good approximation to the wave function of 
the CED state. In sharp contrast to the usual Lanczos method, 
by taking F around the energy of this state, the filter diagonal- 
ization handles only the target CED state. This is the reason 
why the filter diagonalization can succeed in the calculation. 

CED have been studied in the shell-model framework by 
using the formalism introduced by Zuker et al. Ii25il . In 
this description, the Coulomb interaction is separated into a 
monopole term Vc,„ and a multipole term Vcm- While Vcm 
accounts for single-particle and bulk effects, Vqm contains all 
the rest. The monopole term Vcm is further divided into the 
single-particle correction £//, the radial term Vc, , and the elec- 
tromagnetic spin-orbit term C/^.. It has been reported that the 
radial term Vcr does not give a correct description in the lower 
part of the upper /p-shell fl^, and furthermore, this term may 
be important only for high-spin states [ 1 8.1 . Therefore, Vcr is 
not considered in the present calculation. The isospin noncon- 
serving interaction is also neglected for the upper half of the 
//7-shell region because the fi orbit is almost not active. 

The single-particle energy shi ft £/ , takes into account the 
relativistic spin-orbit interaction [26]. This interaction origi- 
nates from the Larmor precession of nucleons in electric fields 
due to their magnetic moments, which, as well known, affects 
the single-particle energy spectrum. £/, can be written as li26ll 
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where niN is the nucleon mass and Vc is the Coulomb poten- 
tial due to the ^^Ni core. The free values of the gyromagnetic 
factors, = 5.586, = 1 for protons and = -3.828, 
= for neutrons, are used. By assuming a uniformly 
charged sphere, g/^ is calculated in the present work by using 
the harmonic-oscillator single-particle wave functions. De- 
pending on proton or neutron orbits, the shift can have oppo- 
site signs. It depends also on the spin-orbit coupling, as for 
instance {I ■ s) =1/2 when j = I + s and {I -s) = — (/ + l)/2 
when j = I — s. As one will see, these strongly affect the re- 
sults in the following discussion. 

With inclusion of the Vcm, £ii, and £/, terms, shell-model 
calculations using the filter diagonalization method are car- 
ried out in the pf^ /2g9 ji model space for odd-odd N = Z nu- 
clei and their even-even IAS partners for A =66, 70, 74, and 
78. In Fig. |2] the calculated CED are compared with available 
experimental data. As can be seen, the calculation reproduces 
the experimental CED remarkably well. In particular, the neg- 
ative CED for A = 70 and the large, positive CED for A = 74 
are correctly obtained. 

We now analyze the results by looking at the components. 
Figure [3] shows the total CED denoted by VcM+is+ih ™d the 
separated multipole, spin-orbit, and orbital parts by Vcm, ^Is, 
and £//, respectively. For the A = 66 pair ^^As/^^Ge, the spin- 
orbit component e/^ has a negative value for the 7 = 0, 2 and 
4 states, while the other two, Vcm and £//, are positive. Since 
the positive values cancel out the negative ones, the net CED 
are small and positive. At 7 = 6, the sudden change from this 
behavior can be understood as the fact that the involving states 
are not the first, but the second 6+ states in the A = 66 pair 
The values for the A = 70 pair ^"Br/^'^Se are large and nega- 
tive for £/,, but positive for Vcm and nearly zero for £;/. Since 
the absolute values of £/, are larger, the total CED are there- 
fore negative. This suggests that the spin-orbit contribution is 
responsible for the observed negative CED in '"^BrP'^Se. For 
the A = 74 pair ^^Rb/^^Kr, the components indicate a similar 
overall behavior as those of the lower spin states in the A = 66 
pair. However, both Vcm and £// are found larger in the A = 74 
pair, and thus are dominant in the summation. Therefore, the 
A = 74 CED are large and positive. For ^^Y/^^Sr, all com- 
ponents are small, and the total CED therefore indicate small 
and positive values. 

Our calculation suggests that for all the pairs studied, the 
spin-orbit component, £/j, is negative. The effect is particu- 
larly enhanced for A = 70. Thus the spin-orbit term appears 
to be the origin of the observed CED anomaly in ^"Br/^^Se. 
But why the absolute values of the spin-orbit component are 
anomalously large only for ^"Br/^^Se? We previously men- 
tioned that the spin-orbit term affects single-particle correc- 
tions differently for neutrons and protons. We may discuss 
these with the following estimate. Due to this interaction, the 
proton gg 12 orbit is lowered by about 66 ke V, while the or- 
bit is raised by about 66 keV. Then the energy gap between the 
proton ggj2 and /5/2 orbit decreases roughly by 132 keV. On 
the contrary, the neutron gg /2 orbit is raised by about 55 keV 
and the orbit is lowered by about 55 ke V, and therefore. 
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FIG. 2: (Color online) Comparison of calculated CED with experi- 
mental data for mass number A =66, 70, 74, and 78. Note that for 
A =66, the calculated CED correspond to the second 6+ states. 



the spin-orbit contribution enhances the neutron shell gap be- 
tween the g9/2 and f^j^ orbit roughly by 110 keV. The total 
effect thus amplifies the difference between the neutron and 
proton orbits by about 242 keV. 

By introducing the difference in neutron occupation num- 
ber between the ground state of Oj'^ and the excited state with 
spin J 



i'9/2 



S'9/2' 



(6) 



it is straightforward to see that the excited states in ^"Se lie 
higher than those in ^"Br with the same spin, and therefore, 
CED in Eq. ([T]l become negative. Since the energy gap be- 
tween the ^9/2 and fp shell for neutrons is by 242 keV larger 
than that for protons, the spin-orbit components can be esti- 
mated to be -33.4, -58.8, and -77.2 keV for spin 7 = 2, 4, and 




FIG. 3: (Color online) Decomposition of theoretical CED for states 
shown in Fig. |2]into three terms (see text for explanation). The total 
CED are denoted by black square symbols. 
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FIG. 4: (Color online) Difference in the neutron occupation number 
,^ between the ground state and the excited states (defined in Eq. 



1^) with spin / =2, 4, and 6 for ^^Ge, ™Se, ^^Rb, and ^'^Y. 



6 from (An^^^^ (™Br) - An^^^^ (™Se)) x 242 keV, respectively. 
These values agree well with the calculated spin-orbit compo- 
nents for ™Br/™Se in Fig. [3] 

Figure|4]shows the differences in neutron occupation num- 
ber, Ani , , between the ground 0,^ state and the excited J- 

' Sg/2 1 

States for A = 62, 66, 70, and 74. One can see that the 
numbers are larger for A = 70. The enhanced neutron exci- 
tations from the fp orbits into the /2 orbit combined with 
the single-particle energy shifts C/^ due to the spin-orbit inter- 
action are the major factors to explain the large negative CED 
in ^*'Br/^*'Se. At the mass number A = 70, the increase of the 
^9/2 occupation reflects the enhanced excitations from the fp 
shell into the gg orbit, which means that ^'^Br and ^'^Se are 
transitional nuclei. In fact, it has been known that the N = Z 
nuclei around A — 70 such as ^^Se and ^^Kr do not correspond 
to a single well-developed shape, and may exhibit the phe- 
nomenon of oblate-prolate shape coexistence 

lliOl- In addi- 
tion, a shape phase transition with an abrupt change in struc- 
ture when the proton and neutron numbers cross N = Z = 35 
has been suggested ff\. 

For ^"Se, our calculation with the standard effective charges 
1 .5e for protons and 0.5e for neutrons obtains B{E2,I ^ ^ ~ 2) 
values of 345, 539, and 594 e^fjn^ for E2 transitions de- 



exciting the 



2+ 4+ 



and 6j states, respectively. These 



are in a good agreement with the experimental values 342 
(19), 370 (24), and 530 (96) e^fm'^ found in Ref. HI. Our 
calculated spectroscopic quadrupole moments are 37.5, 50.2, 
and 55.5 efm^, respectively for the — 2\,A^^ and 6j'^ states. 
Using these values and assuming an axial deformation, we 
may estimate the quadrupole deformation as —0.21, —0.22, 



and -0.22 for the r 



2+ 4+ 



and 6[ states, respectively. 



suggesting that the ^''Se yrast states are oblately deformed. 
This is in consistence with the conclusions from the recent 



measurement 11151] . The deformation is roughly constant, and 
does not show changes with increasing spin. Thus from our 
shell-model calculation, the negative CED for A = 70 are 
not attributed to subtle differences in the Coulomb energy as 



shapes evolve with spin Jstl ■ 

In conclusion, by performing modern shell-model calcu- 
lations, we have investigated the CED effects between the 
isospin r = 1 states in odd-odd N — Z nuclei and the iso- 
baric analogue states in their even-even neighbors for the up- 
per /p-shell nuclei. In order to obtain the T = \ states for 
odd-odd N = Z nuclei, we have gone beyond the usual Lanc- 
zos method by employing the filter diagonalization method 
for the first time in application. It has been shown that the 
anomalous CED found for the pair ^"Br and ^"Se originates 
from neutron excitations from the //?-shell to the g<) /2 intruder 
orbit reflected in a sudden enhancement in the electromag- 
netic spin-orbit term for A — 70. The study has shown how 
the structure changes manifest themselves in a measure of the 
isospin-symmetry breaking in effective nuclear interactions. 

This work was partially supported by the National Natural 
Science Foundation of China under contract No. 11075103 
and 11135005. 
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